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The straightforward method of transformation optics implies that one starts from the coordinate 
transformation, determines the Jacobian matrix, the fields and material parameters of the cloak. 
However, the coordinate transformation appears as an optional function: it is not necessary to 
know it. We offer the solution of some sort of the inverse problem: starting from the fields in the 
invisibility cloak we directly derive the permittivity and permeability tensors of the cloaking shell. 
The approach can be useful for finding material parameters for the specified electromagnetic fields 
in the cloaking shell without knowing coordinate transformation. 

PACS numbers: 81.05.Xj, 78.67.Pt, 41.20.Jb 
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I. INTRODUCTION 



Precursor theoretical papers [1HJ] and the rapid devel- 
opment of the metamaterials as technological branch of 
optics and material physics had led to the discovery of 
the exciting transformation optics Q . The most interest- 
ing application of the transformation optics is obviously 
the invisibility cloaking. Some groups recently reported 
on the realization of the experimental invisibility of large 
objects @, 0| using the carpet cloaking [1, Transfor- 
mation optics approach is extended now to describe the 
field concentrators perfect lensing [ll|, illusion op- 
tics fl2rll3 | , optics in curved space [lj| , optical analogues 
of the cosmological redshift [la ] , Aharonov-Bohm effect 



and black hole horizon [11 1 



In general, transformation optics is an approach deal- 
ing with the change of material parameters and field dis- 
tributions owing to the transformation of a spatial region. 
The transformation is usually a sort of compression or 
tension. It can be described by the coordinate transfor- 
mation converting electromagnetic (original) space to the 
physical (transformed) space. Since the objects to hide 
are commonly in the free space, the original space is the 
vacuum. The squeeze of vacuum of electromagnetic space 
corresponds to the shell in physical space. Spatial region 
inside the shell does not interact with the incident radi- 
ation. Therefore, it is not visible and can be filled with 
any material or contain any object inside. The cloaking 
shell is characterized by the very complex material pa- 
rameters: (i) both dielectric permittivity and magnetic 
permeability tensors should be equal, and (ii) material 
should be anisotropic and inhomogeneous. Nevertheless 
the invisibility cloaking has been already experimentally 
verified |l[l7Hli. 

There are several methods to derive the expressions 
of transformation optics. Except the classical [f| and 
general-relativity-based [ll[ approaches the alternative 



methods can be mentioned. Primarily it is the approach 
of Tretyakov et al [2l[ , which introduces the general con- 
cept of field transformations defined as linear relation be- 
tween the original and transformed fields. The transfor- 
mation optics of Yaghjian et al [22j uses the field bound- 
ary conditions, but not the boundary conditions embed- 
ded into the coordinate transformation as usual. The 
so called method of generating functions (or inverse ap- 
proach) [H, [24| is offered for some special cases of spher- 
ical and cylindrical cloaks. The generating function (for 
instance, some component of the permittivity tensor) is 
used instead of the coordinate transformation. 

The approach of this paper is not intended for the re- 
placement of the ordinary method of transformation op- 
tics. The aim of the approach is to find the cloaking 
parameters of the material for the predefined fields. The 
reconstruction of the medium parameters from the field 
distribution can be called an inverse problem in trans- 
formation optics by analogy with the inverse problem in 
scattering theory. On the one hand, we use the known 
results of transformation optics, such as the formulae for 
the fields and material parameters expressed in terms of 
the Jacobian matrix. On the other hand, the coordinate 
transformation in explicit form is not necessary. In the 
approach, we specify the fields inside the cloak assuming 
the incident plane wave. Then the dielectric permittivity 
and magnetic permeability tensors are derived using the 
three scalar field potentials describing the wave phase, 
electric field vector, and magnetic field vector. 



II. INVERSE PROBLEM 

Let the electric and magnetic fields in the electromag- 
netic and physical spaces are specified. We are look- 
ing for the Jacobian matrix J and material parameters 
e and fi corresponding to such fields. Electromagnetic 
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FIG. 1: Sketch demonstrating the geometry of the problem. 



fields in electromagnetic E'(r'), H'(r') and physical E(r), 
H(r) spaces are assumed to be functions of the coordi- 
nates of their own spaces, that is r' = (x[, x' 2 , x' 3 ) and 
r = {x\, x<ii X3), respectively. In general, the Jacobian 
matrix is the function of coordinates in both electromag- 
netic and physical spaces, because J(r, r') connects the 
coordinates r and r'. The link between the coordinates 
(coordinate transformation) can be found using the equa- 
tion dxj/dx'i = Jij. Then the Jacobian matrix can be 
written as a quantity of the single space, J(r) or J(r'). In 
this paper we are not interested in the coordinate trans- 
formation. Below we will make sure that the coordinate 
transformation is not necessary for finding dielectric per- 
mittivity and magnetic permeability of the transformed 
medium. 

In electromagnetic space, the plane wave propagates 
in z'-direction in vacuum e' = 1 and \i' = 1. Fields 
of the plane wave are described in the common way: 
E' = exp(ikoz')e x an d H' = exp(ifcoz') e y ( see figure [TJ, 
where fco = w/c is the wavenumber in vacuum, oj is the 
circular frequency, and c is the speed of light in vacuum. 
So, we consider the transformation from the plane wave 
in vacuum to the complex field distribution in physical 
space. Since an arbitrary wave can be presented as the 
superposition of plane waves, we expect that the invisi- 
bility cloaking for the plane wave will give the material 
parameters applicable for any field distribution in the 
electromagnetic space. 

Boundary conditions are usually specified via the coor- 
dinate transformation: the spatial region in electromag- 
netic space confined by the surface Si is squeezed to the 
region between the surfaces Si and S 2 - If the coordi- 
nate transformation is not defined, the boundaries of the 
transformed space should be set by means of the fields. 
We have four "no-scattering" conditions (the continuity 
of the electric and magnetic fields tangential to the outer 
interface Si defining the impedance matching with the 
ambient medium) and the "no-field-inside-cloak" condi- 
tion, which can be formulated as zero energy flux passing 



through the inner boundary (field does not penetrate): 

x E(Si) = ni X E'(Si), m x H(Si) = n x X H'(Si), 
n 2 S(S 2 ) = -^n 2 [E(S 2 ) x H(S 2 )] - 0, (1) 

47T 

where ni and n 2 are the unit vectors normal to the sur- 
faces S\ and 5*2, respectively. For example, the interfaces 
S\ and 5*2 can be spheres or cylinders. But more com- 
plicated geometries are possible, too. We emphasize that 
only the inner and outer interfaces are specified by equa- 
tion (|T|), but not the coordinate transformation. 

Further we will find the inverse Jacobian matrix from 
equations of transformation optics 

E(r) = J-\r, r')E'(r'), H(r) = J-\r, r')H'(r'). 

(2) 

Solution of the first vector equation can be written as 
the sum of the contributions along the electric field E' 
and orthogonal to the electric field E': 



J- 



E® E' 
E' 2 



(3) 



where a is an arbitrary matrix. Dyad E ® E' (in in- 
dex form (E ® E')jj = EiE'j) operates as a projec- 
tor for arbitrary vector a: (E ® E')a = E(E'a) and 
a(E <E> E') = (aE)E'. Three-dimensional antisymmetric 
tensor E /x dual to the vector E' (25l . I2H ] can be repre- 
sented in index form as (E' x )ife = J2^=i £ijkEj, where 
Eijk is the Levi-Civita symbol. Dual tensor is charac- 
terized by the following properties: E' x a = E' x a and 
aE' x = a x E'. By substituting solution (|3]) into the sec- 
ond equation in ^ we obtain the equation with respect 
to matrix a: 



H 



(WW) 
E' 2 



E = d(E' x H'). 



(4) 



The solution of this equation can be presented in the 
similar form as ©. Using the relationship P' X E' X = 
E' (g> P' we finally derive the Jacobian matrix 



J- X (r, 



E <g> H' - H <g> E' 

p/2 



P ,x +a®P', (5) 



where quantity P' = E' x H' is proportional to the Poynt- 
ing vector S' of the plane wave in electromagnetic space. 
An arbitrary vector a (three elements of matrix J) can- 
not be found from equation @. 

Actually we have found more than we expected, be- 
cause from expression © it follows that the Jacobian 
matrix remains the same, if the fields in electromagnetic 
and physical spaces are simultaneously multiplied by the 
same scalar function f'(r') = f(r). This means that the 
Jacobian matrix is invariant with respect to the phase 
factor. If E' = exp(ikoz')e x and H' = exp(ifcoz') e y) then 
the fields in physical space are E = exp(ifcoz'(r)) J~ 1 e x 
and H = exp(ikoz' (r))J~ 1 e y . Therefore, we can spec- 
ify just fields E' = e x and H' = e y ignoring the same 
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scalar functions, which are terminated in the definition 
of Jacobian matrix. So, we derive 

J- 1 (r) = E ®E{, + Hi)®H^ + ao®Po > (6) 

where E = e x , H = e y , P' = e z , E = exp(-ifc z')E, 
Hq = exp(— ifcoz')H, and ao = exp(2ifco^')a. 

Vectors E , H , and a of the physical space are not 
arbitrary. In addition to the boundary conditions ([T]) 
restricting the electric and magnetic fields, the Jacobian 
matrix J is not arbitrary as well. It should be agreed with 
Maxwell's equations or, in terms of the transformation 
optics, the identity condition 

VxJ- x (r) = (7) 

should hold. Introducing matrix ([5]) to equation ([7]) we 
get to tree vector relationships 

V x E = 0, V x H = 0, Vxaj = 0. (8) 

Potential vectors Eo, Ho, and ao defined by ((SJ) can be 
presented by means of the scalar potentials ip, Xt an d V 

as 

E - Vip, Ho = Vx, a = V77, (9) 

respectively. The meaning of the potentials can be de- 
termined from the comparison of the Jacobian matrix 

J 1 = V <8> (ip&x + X e y + V e z) (10) 

with its definition J -1 = V (S> r'. The correspondence 
is clear: scalar potentials are the coordinates in elec- 
tromagnetic space, that is x' — ip(r), y' — x( r )) an d 
z' = T)(r). These three equations define the coordinate 
transformation. Since the plane wave in electromagnetic 
space propagates in the z'-direction, the phase factor is 
described by the potential r\ as exp(ifcoz') = exp(ifcor/). 
The potentials have clear physical meaning for the fields: 
ip and x define the vector structure of the electric and 
magnetic fields, while r\ is the phase function. Thus, the 
electric and magnetic fields in the cloaking shell are 

E = e ifeo ' 7 WV0(r), H = e ifcn " (r) Vx(r). (11) 

where ip, Xj an( i V are the electric, magnetic, and phase 
potentials, respectively. 

Dielectric permittivity and magnetic permeability ten- 
sors equal 

A(r) = Vip ® Vip + Vx ® V% + V77 ® Vrj, (12) 

where the superscript T denotes the transpose operation. 
It is important that the permittivity and permeability 
tensors are defined in the physical space. That is why the 
coordinate transformation is not required and the tensors 
are set in terms of the potentials — field characteristics. 



It should be mentioned that we do not differ the covariant 
and contravariant indices in three-dimensional tensors. 
We use the tensors in index- free form, e.g. e — (ey) = 
(e 13 ). The methods of operating with index- free objects 
are well developed in Refs. [25l - l27| and applied in the 
optics of complex (anisotropic, bianisotropic) media [28l — 

m 

One more limitation on the potentials is caused by 
the fact that the Jacobian matrix is nonsingular, i.e. 
det(J _1 ) ^ or det(A) ^ 0. In other words, the vec- 
tors Eo, Ho, and ao are not parallel: (Eo x Ho)ao / 0. 
Since Eo and Ho are not parallel, vector ao should surely 
have the component along Po = Eo x Ho. 

In general, all potentials are important for setting elec- 
tric and magnetic fields in the cloaking shell. However, if 
only vectorial distributions of the electric and magnetic 
fields are essential for us, the phase of the field can be 
omitted. Then the same vectorial distributions follow for 
various phase potentials and we can tune the phase po- 
tential to find the simplest form of the permittivity and 
permeability tensors. 

How should the potentials look like? We can sup- 
pose that they should be like Cartesian coordinates writ- 
ten in curvilinear coordinates, because the meaning of 
the potentials are the Cartesian coordinates in the elec- 
tromagnetic space. For example, potentials ip = x' = 
R(r) sin 9 cos ip, \ — y' — R(r) sin 8 sin ip, and T) = z' = 
R(r) cos 9 can be applied for determining material ten- 
sors in spherical coordinates. Detailed discussion of the 
possible forms of the potentials is given in the subsequent 
sections. 



III. SPHERICAL CLOAK 

Let us consider a spherical cloak with inner and outer 
boundaries a and 6, respectively. External electric and 
magnetic fields (they coincide with the fields in the elec- 
tromagnetic space) 

E' (r) = e lfc ° r cos 9 (sin 9 cos ipe r + cos 9 cos (peg - sin ipe^ ) , 
H'(r) = e lfc ° rcose (sin 9 sin tpe r + cos9smipe e + cos^e^) 

are linked with electric and magnetic fields inside the 
cloak 

E(r)=e^(-) f^ + ige. + ^lV), 
\or rod rsm9otp J 

HW = e^) de r + i|e e + ^-^(13) 
\or r 08 r sin dip J 

via the conditions of continuity of the tangential fields 
(HJ) at the outer boundary r — b. These conditions can 
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be written in the form 



d^{b,9,p>) 

cos 9 — r)(b,0, </?)) ^ 



COS 



09 

dx(b,0,ip) _ jk(b cos e-r,(b,e 



dip 



,v)) 



sin v cos (p. 



(14) 



For the agreement of the equations a couple of par- 
tial differential equations d 2 ip/d9dp = d 2 ip/dpd9 and 
d 2 x/dOdip — d 2 x/dipd9 with respect to the func- 
tion 77(6, 9, ip) should hold true. In the current sit- 
uation, this couple of equations is easily solved giv- 
ing rj(b,9,p) = bcos9 + 770. Electric and mag- 
netic potentials at r = b follow from the integration: 
ip(b, 9, p) — 6exp(— ifc?7o) sin 9 cos p + ipo and \(b, 9, ip) — 
6exp(— i/c?7o) sin6*siniy9 + xo, where 770, ipo, and xo are 
constants. It can be an arbitrary dependence of the po- 
tential on the radial coordinate r reducing to ip(b,0,tp), 
x(b, 9, ip), and 77(6, 9, ip) at r — b. From the great variety 
of functions we choose 
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FIG. 2: Differential cross-sections for spherical cloaks 
with different field profiles /(?"): linear f(r) = b(r — 
a)/(6 — a), quadratic f(r) = b(r 2 
soidal f(r) = 6sin[(7r/2)(r - a)/(6 
fir) = (26/tt) arctan[(r - a)/(b - r)] 
kob = 2-7T. 



- a )/(& — a ), sinu- 
a)], and arctan-form 
Parameters: koa = n, 



ip(r,9,p) = f(r) sin 9 cost/?, x(r, 9, ip) = f/(r) sin 9 simp, 
r)(r,9,<p) = h(r)cos9, (15) 

where /, g, and h are arbitrary functions except ful- 
filment of the following boundary conditions: f(b) — 
g{b) = h(b) = b, and f(a) = or g(a) = 0. Two 
last conditions follow from the vanishing of the Poynt- 
ing vector component normal to the inner sphere r = a: 
S r ~ /( a ).9( a ) — 0. The cloak is realized for any /, g, 
and h discussed above. So, we can choose these func- 
tions to provide the needed electric and magnetic field 
vectors (and phase, if required). Then the permittivity 
and permeability tensors are easily calculated from equa- 
tion (fHj) . 

If the phase is not important, potential 77 can be used 
for optimization of the permittivity and permeability ten- 
sors. Assuming f(r) — g(r) ^ h(r) we derive 

fP fP T 

£ = ~^f e r ® e r + —e e <%>e g + -je v ® e v 

, f cos 6(h'h-f'f) , 

H — fe r eg) e v + e v ® e r ) , (16) 

ri 

where P(r, 6») = / 2 cos 2 6»+ft 2 sin 2 (9, T(r, 9) = h'f cos 2 (9+ 
f'h sin 2 6*, /' = df/dr. The same distribution of the 
electric |E| 2 and magnetic |H| 2 fields and time- averaged 
Poynting vector (c/87r)Re(E x H*) is realized for any 
phase function h(r). Generally dielectric permittivity 
tensor depends on two coordinates, r and 9, and has non- 
diagonal elements. The simple phase functions ft. = 6 or 
h = r do not eliminate the ^-dependence or anisotropy. 
They can be eliminated only for the equal radial func- 
tions, i.e. f(r) — g(r) — h(r). Then the permittivity 
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FIG. 3: Magnetic field component of the arctan-form 
cloak for <p = and 9 — 45° calculated both analytically 
and using the simulation. In the insets the magnetic field 
distributions in the plane ip = are shown. Parameters: 
koa = 7r, kob = 2-7T. 



tensor takes the well-known form 
f 2 

e = ~^rp e r ® e r + / eg <S> ee + f e v <g> e^. (17) 

Example of the more sophisticated choice of the po- 
tentials is 

tp{r, 9, p) = f(r) sin 0(r, 9) cos $(r, p) + ipo{r), 
X(r,9,p) = g(r) sin6(r, 0)sin$(r, ip) +Xo[r), 

rf(r,e,(p) = ft(r) cos e(r,0), (18) 
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FIG. 4: Magnetic fields H y evaluated (a) analytically and 
(b) in simulation. Electric field component E x obtained (c) 
analytically and (d) using the simulation, (e) Magnetic H y 
and (f) electric E z fields simulated in the case of the plane 
wave incident in x direction. Parameters: koa = tt, kob = 2-n. 



where f(b) = g(b) = h(b) = 6, 9(6, 6) = 9, and $(&, ip) = 
<p. 

Any function f(r) — g(r) — h(r) meeting the appro- 
priate boundary conditions (f(b) = b and f(a) = 0) can 
be used to obtain the cloaking parameters (IT71) . In fig- 
ure [H the differential cross-sections for four field profiles 
f(r) are shown. Inhomogeneous cloaking shell is divided 
into N — 30 homogeneous spherical layers. Differential 
cross-section normalized by the geometrical cross-section 
of the spherical particle is calculated using the analytical 
expressions given in Ref. [3lJ. As it has been expected, 
all curves in figure [5] show good cloaking properties in 
the whole angle range. 

In figure[3J we compare two solutions for the cloak with 
arctan profile: one curve is obtained as the closed-form 
solution with f(r) — (2b /w) arctan[(r — a)/(b — r)], an- 
other curve is computed for the appropriate cloaking pa- 
rameters using the field scattering technique. The fields 
in the cloak are very close one to another. The only re- 
gion, where they differ, is the region inside the cloaking 
shell. It is explained by the inaccuracy of the simulation 
which cannot precisely evaluate the fields near r = a (the 
hidden region is filled with glass). 



IV. CYLINDRICAL CLOAK 

The similar algorithm can be applied for the cylindrical 
cloak extended from r — a to r = b (cylinder axis is y- 
directed). 

(i) Write down the electric and magnetic fields in and 
outside the cloak. Fields in the cloaking shell are given 
by equation pTjl . 

(ii) Apply the boundary conditions ([TJ) at the outer Si 
and inner S2 interfaces. 

(iii) Choose the needed form of the potentials ip, Xi 
and 77 and determine the fields (fTTj) . 

(iv) Calculate the permittivities using equation (|12l) . 
We pick out the potentials for the cylindrical cloak in 

the form 

tp(r, ip, z) = f(r) cos (p, \( r , <P, z ) = 9( r )v, 

r](r, ip, z) = h(r) sin ip. (19) 

Then the fields are equal to 

E = e *>*(r)a« v (|; cos per - ^sin^v) , 

H = c lkoh ^ sin v(^£.ye r + g(r)e y ^j , (20) 

where f(b) — h(b) = b, g(b) = 1 and f(a) = or g(a) = 0. 
The material parameters of the ordinary cylindrical cloak 
can be obtained for f(r) — h(r) and g(r) = 1. 

If / 7^ h, we can expect more complicated material ten- 
sors. Magnetic field in the cloak equals H = exp(ifc?7)e a . 
It possesses the unit amplitude and can differ only in 
phase. We take f(r) = b(r — a)/(b — a) and h(r) = b 
satisfying the required boundary conditions. Function / 
shall be chosen to provide the needed electric field. The 
phase potential in our case is very simple: 7/ = b sin ip. 
Nevertheless, the permittivity tensor is not diagonal in 
the basis of the cylindrical coordinates: 




where e rr = [(b — a) 2 + (r — a) 2 t&nip] / [(b — a)r], e w — 
r/(b — a), E yy — b 2 cos 2 ip]/[(b — a)r], and e rtp — (r — 
a) tan <£>/(& — a). In figure SJa) and (c) we demonstrate 
the fields, which we expect in the cloak. They are plot- 
ted using the closed-form expressions H y = exp(ifc& sin ip) 
and E x = exp(ifc6sin<p)(l — a sin 2 ip/r)b/(b — a), respec- 
tively. To verify that we will obtain the similar field 
distributions in the cloaking shell with parameters (|21l) . 
the numerical simulation was carried out. From figures 
SJb) and (d) it is obvious that the cloaking shell indeed 
guides the prescribed fields. However, it will be another 
field distribution, if the incident wave propagates along 
x axis (or in some direction differing from z axis). In 
figures H] (e) and (f ) exactly such a situation is on view. 
One indicates the cloaking property of the shell in this 
case, too. 
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V. CONCLUDING REMARKS 

The inverse approach discussed in this paper can be 
applied to other transformation optics devices. For ex- 
ample, concentrator consists of the inner and outer re- 
gions [HI H|[ • The inner region follows from the squeeze 
of the electromagnetic space, while the outer one is the 
result of extension. Both transformed spatial regions are 
characterized by the parameters (fTTI) , but with their own 
potentials: tpi, Xi> Vi an d V^, X2, V2 for outer and inner 
regions, respectively. The potentials are not arbitrary. 
They arc restricted by the boundary conditions put on 
the fields. The impedance matching conditions are the 
first two vector equation in ([1} . Two more equations are 
the continuity of the tangential components at the inter- 
face between the inner and outer regions. The field will 
be concentrated in the inner region, if there is a squeeze, 
i.e. condition det(J _1 ) > 1 holds true. 

The second example for application is the cylindrical- 
to-plane- wave convertor [33j • The idea of this device is 
to modify the cylindrically symmetric region to the part 
of a square. Then the cylindrical electromagnetic wave 
generated at the center will turn to the plane wave af- 
ter transmission through the transformed region. In our 
notations, there is a part of square with the plane wave 
electromagnetic field outside. The fields inside the square 
have the form . We put only the impedance- matching 
boundary conditions (JXJ) on the fields to obtain the con- 



vertor. 

In conclusion, we have solved the inverse problem in 
transformation optics, that is we derived the analyti- 
cal formula (| 12[) for the material tensors using the in- 
formation about the fields inside the cloaking shell: the 
phase and vector distributions of the electric and mag- 
netic fields. The closed-form expressions are not gen- 
eral. They are written for the incident plane wave and 
e' = // = 1 in the electromagnetic space. Nevertheless, 
there is no limitations on the form of the cloak. More- 
over, the derived material tensors describe the cloak for 
any incident wave, not only for the plane wave used to 
obtain s and /1. If any other wave is incident, the cloak- 
ing property of the shell keeps, but the field distribution 
in the cloak is not specified now. It is well known that 
material parameters of a wide class of cylindrical cloaks 
have a singularity at the inner boundary. We believe 
that the inverse approach can be useful in finding non- 
singular parameters of the cloaking shell. The example 
of nonsingular parameters is equation (|21[) . The method 
developed seems to be very attractive as the means of 
deeper understanding of the physics and mechanisms of 
the optical cloaking and can substantially help in design 
of the cloaks with predefined fields inside. 

Financial support from the Danish Research Council 
for Technology and Production Sciences via project THz 
COW and Basic Research Foundation of Belarus (grant 
F10M-021) is acknowledged. 
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